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^1 ABSTRACT 

riZ| We perform hydro dynamical simulations of minor- merger induced gas sloshing and the 

&* subsequent formation of cold fronts. Using the Virgo cluster as a test case, we show 

q for the first time that the sloshing scenario can reproduce the radii and the contrasts 

Vh in X-ray brightness, projected temperature, and metallicity across the observed cold 

^ fronts quantitatively. The comparison suggests a third cold front 20 kpc NW of the 

C3 Virgo core. 

We identify several new features typical for sloshing cold fronts: an alternating dis- 
tribution of cool, metal enriched X-ray brightness excess regions and warm brightness 
deficit regions of reduced metallicity; a constant or radially decreasing temperature 
^^ accompanied by a plateau in metallicity inside the cold fronts; a warm rim outside the 

^^ cold fronts; and a large-scale brightness asymmetry. We can trace these new features 

{N| not only in Virgo, but also in other clusters exhibiting sloshing cold fronts. 

By comparing synthetic and real observations, we estimate that the original minor 
■^ merger event took place about 1.5 Gyr ago when a subcluster of 2-4 x 10 13 Mq passed 

^^ the Virgo core at 100 to 400 kpc distance, where a smaller mass corresponds to a 

£^ smaller pericentre distance, and vice versa. From our inferred merger geometry, we 

t-H derive the current location of the disturbing subcluster to be about 1-2 Mpc E of the 

Virgo core. A possible candidate is M60. 

Additionally, we quantify the metal redistribution by sloshing and discuss its im- 
portance. 
;_5 We verify that the subcluster required to produce the observed cold fronts could 

be completely ram pressure stripped before reaching the Virgo centre, and discuss the 
conditions required for this to be achieved. 

Finally, we demonstrate that the bow shock of a fast galaxy passing the Virgo 
cluster at ~ 400 kpc distance also causes sloshing and leads to very similar cold front 
structures. The responsible galaxy would be located about 2 Mpc north of the Virgo 
centre. A possible candidate is M85. 

Key words: galaxies: clusters: individual: Virgo - galaxies: clusters: individual: A496 
- galaxies: clusters: individual: Perseus - galaxies: individual: M87 X-rays: galaxies: 
clusters methods: numerical 



1 INTRODUCTION medium (ICM) of galaxy clusters, among them cold fronts. 

The first ones were noticed as sharp discontinuities in X-ray 
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During the last decade, high-resolution X-ray observations , . , , . , £ ,, r A01/(n jtt — ; r— ;- 

. ° . . iir i ■ i brightness on one side of the cores or A2142 ( Markevitch 

have revealed a wealth of structure in the mtra-cluster 



tetal. 2000), A3667 (Vikhlinin et al. 2001), and the bullet 



cluster IE 0657-56 (Markevitch et al. 2002). Temperature 
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measurements revealed that these brightness edges were not 
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the expected bow shocks but cold fronts (CFs), where the 
brighter and denser side is also the cooler one. For a shock, 
the temperature jump across the brightness edge would be 
opposite. Several such CFs associated with merging clus- 



ters have been found (see review by Markevitch & Vikhlinin 
2007| , hence they are called merger CFs. They have been in- 
terpreted as the contact discontinuity between the gaseous 
atmospheres of two different clusters. More recently, also 
shocks associated with the merging have been discovered in 
some of these clusters (see Markevitch|2010 for a review) . 

Additionally, a second class of CFs, named sloshing 
CFs after their most likely origin, have been detected. Here, 
the CFs form arcs around the cool cores of apparently re- 
laxed clusters. This type of CFs is the subject of this pa- 
per. For clarity, we will first introduce the sloshing scenario 
(Sect.ITjJ) which is commonly believed to explain this type 
of CFs. Section 11.21 summarises the observational charac- 
teristics of sloshing CFs known so far. Targeting the CFs in 
the Virgo cluster, we show that the sloshing scenario can ex- 
plain the observational characteristics not only qualitatively, 
but also quantitatively. Furthermore, we obtain insights in 
the correct way of interpreting CF observations, and we un- 
cover a number of new observational characteristics, which 
can also be transferred to other clusters. 



1.1 Gas sloshing scenario 



Originally proposed by Markevitch et al. (2001) (see also 
review by iMarkevitch & Vikhlinin 2007), the idea of the 



sloshing scenario is the following: A gas- free sub cluster falls 
through the main galaxy cluster, which initially has a near- 
hydrostatic ICM distribution. During the pericentre pas- 
sage, the interaction slightly offsets the gas in the main 
cluster core, but does not disrupt the main cluster core. 
The offset mechanism is described in detail in Ascasibar & 



Markevitch] (|2006) (AM06 hereafter): while the subcluster 
approaches the main cluster, the core of the main cluster, 
i.e. both, dark matter (DM) peak and cool gas peak, move 
towards the subcluster. When the subcluster passes the main 
core, it gravitationally pulls the main cluster core along. 
However, the ICM velocity field surrounding the main clus- 
ter core opposes this pull. While the DM peak can still move 
towards the receding subcluster, the gas peak is held back. 
After the subcluster has passed the central region and moves 
away, the offset ICM falls back towards the main cluster 
DM peak and starts to slosh inside the main potential well. 
Generally, the highest frequencies of sloshing occur at small 
radii. Thus, gas at a given radius moving into one direction 
will frequently encounter gas from larger radii still moving 
in the opposite direction. Such opposing flows lead to the 
formation of density, temperature, and consequently X-ray 
brightness discontinuities. The formation of a similar dis- 
continuity occurs at the leading side of a cool cloud moving 
through a hotter ambient medium ( Heinz et al.|2003 ). As the 
central cool gas moves outwards into regions of lower pres- 
sure, adiabatic expansion enhances the temperature contrast 
across the fronts. Thus, sloshing CFs are contact discontinu- 
ities between gases of different entropy, originating from dif- 
ferent cluster radii. Usually, the subcluster passes the main 
cluster core at some distance, it transfers angular momen- 
tum to the ICM, and the sloshing takes on a spiral-like ap- 
pearance. 



AM06 have performed hydrodynamical SPH+Af-body 
simulations confirming this scenario. Clearly, the gravity of 
the subcluster does not only influence the main cluster's 
ICM, but also its DM distribution. However, the interaction 
with the subcluster induces only a slow motion of the central 
DM peak w.r.t. the overall cluster potential, while in the 
meantime the ICM sloshes inside the slowly moving central 
potential well. These simulations also demonstrated that the 
sloshing reproduces the morphology of observed CFs, where 
less massive subclusters cause weaker signatures. 

AM06 also tested the impact of a subcluster containing 
a gaseous atmosphere of its own. They showed that the gas- 
gas interaction tends to leave clear observable signatures like 
a tail of ram-pressure stripped gas and an overall much more 
disturbed appearance. As the clusters with observed slosh- 
ing CFs do not show such signatures, the scenario favours 
gas- free subclusters. | Johnson et aL| |2010) presented obser- 
vations of the double cluster Abell 1644, where the two sub- 
clusters A1644N and A1644S appear to be merging. While 
both subclusters contain an ICM of their own, the core 
of A1644S exhibits a characteristic spiral-shaped CF. How- 
ever, in contrast to the clusters with relaxed appearance and 
sloshing CFs, both subclusters of A1644 show a disturbed 
large-scale morphology. 



The simulations of ZuHone et al. (2010) are similar to 



the ones of AM06, but concentrate on the heating efficiency 
of the sloshing process. 

As an alternative to the scenario described above, the 
oscillation of the central galaxy or central dark matter peak 
of the cluster has been proposed (Lufkin et al.|1995 Fabian 



et al.|2001| |Tittley fe Henriksen||2005| . The moving central 
galaxy could experience ram pressure stripping by the sur- 
rounding ICM, and the contact discontinuity between the 
cooler central and warmer ambient gas could form the CFs. 
However, it is unclear whether this scenario would produce 
multiple fronts, which arise naturally in the sloshing sce- 



Birnboim et al. (2010) propose that at least some CFs 



are caused by merging of shocks. In this case, the CFs should 
form quasi-spherical rings around the cluster centre, and 
could be found at large radii out to 1 Mpc. We show that 
sloshing CFs are accompanied by a characteristic large-scale 
asymmetry, which would be absent in this scenario. 

1.2 Observed characteristics of sloshing CFs 

Most sloshing CFs manifest themselves as brightness edges 
forming arcs around the cool cores of clusters that otherwise 
appear relaxed and show no obvious signatures of recent 
merging ( |Markevitch fe Vikhlinm||2007| [Owers et al.||2009 ). 
Often, pairs or even triplets of edges appear on more or less 
opposite sides of the cluster core. Their morphology depends 
mainly on the angle between our line-of-sight (LOS) and 
the orbital plane of the subcluster: If the interaction is seen 
face-on, the edges form a spiral-like structure. If the LOS is 
parallel to the orbital plane, we see arcs on alternating sides 
of the cluster core. The brightness and temperature contrast 
across the edges is more modest compared to the merging 
CFs. Typical temperature contrasts range between 1.2 and 
2.5. 



This type of CF is reported to be ubiquitous (Marke- 



|vitch et al.1|2003l |Markevitch fe Vikhlin in 2007, Ghizzardi 
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et al. 2010), but high-resolution observations are available 



for only 13 clusters: RX J172 0.1+2638 (jMazzotta et al 
2001] |Mazzotta fe Giacintucci| [20081 [Owers| |2008| |Owers 
et al.|2009|), MS145 5.0+223 2 JMazzotta fe Giacintucci|2008| 



Owers et al. 2009), A2142 (Markevitch et al. 2000 lOwers 



A1795 (Markevitch et al. 



Perseus (Churazov et al 



et aI.|2009^ r A496 (|Tanaka et al.||2006| |Dupke let aL|2007 >,' 



A2029 (Clarke et al. 2004 Million & Allen 2009), A2204 



2A0335 +096 (jMazzotta et al.||2003| |Sanders et al ||2009a| ) , 



([Sanders et al.|2005[|2009b[, Ophiuchus ([Million et al.|2009 ) , 



2001 Bourdin & Mazzotta 2008), 



2003 Sanders et al. 2005), Cen 



taurus (|Fabian et al.|2005| |Sanders fe Fabian||2006| ), A1644 



(Johnson et al. 2010). Our work will concentrate on the 



Virgo cluster (Simionescu et al.|2010[ S10 hereafter). 



In addition to the edges in brightness and tem- 
perature, temperature maps show spiral- or arc-shaped 
cool regions inside the actual fronts. In X-ray bright- 
ness residual maps, an excess corresponding spatially 
to the structure in the temperature map is found 
(A2029, A2204, 2A0335+096, Perseus). For some clus- 
ters (RX J1720.1+2638, MS1455.0+2232, A496, Ophiuchus, 
2A0335+096, A2204, Centaurus), even a corresponding 
structure in metallicity maps is seen. As X-ray brightness 
maps can be obtained more easily than temperature or 
metallicity maps, a spiral-shaped excess wrapped around 
the cluster core is often regarded as an indicator for ongoing 
sloshing. 

Figure [I] shows two versions of a residual X-ray bright- 
ness map for Virgo: The top panel displays the brightness 
residuals w.r.t. the azimuthally averaged brightness. The 
bottom panel shows residuals w.r.t. the best-fitting radially 
symmetric /3-model. This map is the same as Fig. 2 in S10, 
but in a colour-scale used consistently throughout this pa- 
per. Both residual maps show a sharp brightness edge about 
90 kpc NW of the core. The edge about 33 kpc towards the 
SE is somewhat more subtle. The Virgo observations are de- 
scribed in more detail in Sect. [3] along with the comparison 
to the simulation results. 



1.3 Aims and outline 

Previous work has shown that the sloshing scenario 
(Sect. [hT]) can in principle explain the known observed char- 
acteristics (Sect. 1.2) of the resulting CFs. The quantitative 
characteristics of the CFs in a cluster (radius, contrast of 
X-ray brightness, temperature, metallicity across the edges, 
morphology in maps) depend on several parameters: 



• host cluster potential and ICM distribution; 

• characteristics of subcluster: mass, size, pericentre of 
orbit, velocity along orbit; 

• age: time since the pericentre passage of the subcluster; 

• angle of our line-of-sight (LOS) w.r.t. the orbit and or- 
bital plane of the subcluster. 

In this work, we concentrate on the Virgo cluster and 
focus on the following aims: 

• Test the sloshing scenario quantitatively. We show that 
we can find a reasonable set of parameters to reproduce 
the observations well (Sect. pi. We note that there is a 
degeneracy between subcluster mass and pericentre dis- 
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Figure 1. Observed residual maps (data divided by average 
or model) of the Virgo cluster. The top panel shows residuals 
w.r.t. the azimuthally averaged brightness. The bottom panel 
shows residuals w.r.t. a /3-profile fit to the observed X-ray bright- 
ness (same as in |Simionescu et aL]|2010| but with colour scale 
consistent to this paper). 



tance/ velocity along orbit, so that only a combination of 
these can be constrained. 

• Refine our understanding of observable signatures of 
sloshing CFs (Sect. [I). 

• Study the ability of sloshing to redistribute metals 
throughout the cluster core (Sect. [5}. 

• Test if the basic assumption of the scenario, i.e., that 
the subcluster is gas- free (see Sect. |1.1[ ), can be true 
(Sect. [6]). Once we have identified a suitable subcluster, 
we show that a gaseous atmosphere could be ram-pressure 
stripped off this subcluster before it reaches the cluster core. 
However, the outcome of this test is not unique, and depends 
on the details of the subcluster's atmosphere. 



• Test an alternative scenario: Churazov et al. (2003) pro- 
pose that the passage of a shock through a cluster can cause 
the formation of CFs and demonstrated this with an ide- 
alised plane shock. We study the more realistic situation of 
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the bow shock of a galaxy and its ability to trigger sloshing 
and CF formation (Sect. [7]). 

To this end, we run hydrodynamical simulations of an 
idealised minor merger of a (gas-free) subcluster with a main 
cluster (see Sect. [2]), which induces sloshing of the ICM of 
the main cluster. We perform a parameter study, varying the 
discussed parameters, and present a detailed comparison to 
the observations (Sects. [3] El. 



2 METHOD 

In a set of idealised merger simulations, we model the in- 
teraction between a massive galaxy cluster ("main cluster") 
and a smaller subcluster ("disturber"). We concentrate on 
minor mergers, where the disturber is significantly less mas- 
sive than the main cluster. A major merger would cause 
additional structure, e.g., destroy the cool core, which is not 
the case for sloshing CF clusters. In the Virgo cluster, only 
in the inner few kpc the cooling time is shorter than 2 Gyr. 
Hence, in the outer regions of interest for our analysis, the 
cooling time is long and we neglect radiative cooling in our 
simulations. Finally, we will only model the first core passage 
of the subcluster and the subsequent gas sloshing. 



2.1 Code 

Our simulations are run using the FLASH code (version 3.2, 
Dub ey et al.||2009 ). FLASH is a modular block-structured 
AMR code, parallelised using the Message Passing Interface 
(MPI) library. It solves the Riemann problem on a Cartesian 
grid using the Piecewise-Parabolic Method (PPM). The sim- 
ulations are performed in 3D and all boundaries are reflect- 
ing. We use a simulation grid of size 3 x 3.5 x 3 Mpc 3 . In most 
simulations, we resolve the inner 32 kpc with Ax — 1 kpc, 
the inner 64 kpc with Ax = 2 kpc, and the inner 128 kpc 
with Ax = 4 kpc. At this moderate resolution, one simula- 
tion can be completed in about 100 CPUh. In Appendix [A] 
we present a high-resolution run of our fiducial case which 
resolves the inner 8 kpc with Ax = 0.5 kpc and the inner 
128 kpc with Ax = 1 kpc, and a very high-resolution run 
with an again twice as good resolution. None of our results 
depends on resolution. 



subcluster. The simplest approximation to account for this 
is to assume that the main cluster responds to the grav- 
ity of the subcluster as a whole, like a rigid body. There- 
fore, we calculate the inert ial acceleration felt by the main 
cluster centre due to the subcluster and add this pseudo- 
acceleration to the entire ICM. This procedure was also 
used by ZuHone et al. (2010). It is a reasonable approx- 



imation for the central region of the main cluster, but it 
is inaccurate for its outer parts, where it will lead to un- 
realistic flows. For small pericentre distances, during peri- 
centre passage of the subcluster the inertial acceleration 
can be large and even produce supersonic motions in the 
outer parts of the main cluster. Hence, we apply the pseudo- 
acceleration only to the central region of the main clus- 
ter inside a characteristic radius, i?dam P , and damp it out- 
side this radius exponentially with a length scale, Ld amp - 
Roediger & ZuHone (in prep.) perform a direct compar- 
ison between sloshing simulations using this rigid poten- 
tial approximation and full hydro+Nbody simulations for 
a more massive cluster resembling A 2029. They show that 
for it^damp > 2a su b, where a su b is the scale radius of the 
subcluster, and Ldamp ~ 0.5i?dam P the rigid potential sim- 
ulations reproduce the full hydro+Nbody simulations well. 
The orientation of the cold spiral structure is reproduced 
accurately. The evolution of the size of the CF structure 
lags behind in the rigid potential approximation by about 
200 Myr. When this lag is taken into account, the size of 
the CF structure is reproduced well. Due to this artefact, 
we will over-estimate the age of the CFs derived from the 
rigid potential simulations by about 200 Myr. For a given 
size of the CF structure, the density contrast across the CFs 
is reproduced well. The asymmetry in density out to at least 
500 kpc is reproduced qualitatively. The temperature at the 
inner side of each CF is also reproduced well. The tempera- 
ture just outside the CFs is too high due to enhanced com- 
pressional heating from unrealistic gas motions outside the 
central region. The damping of the inertial frame correction 
explained above can reduce the over-estimate of the tem- 
perature on the outer side of the CFs somewhat, but cannot 
prevent this artefact completely. 

In the case of Virgo, we use i?dam P = 500 kpc and 
Ldamp = 300 kpc. Apart from lowering the temperature out- 
side the CFs, this choice does not influence the resulting CF 
structure. 



2.2 Rigid potential approximation 

The gas dynamics is described by the hydrodynamical equa- 
tions. Additionally, the ICM is subject to the gravitational 
acceleration due to the main cluster and the subcluster. 
In order to sample the parameter space of the subcluster 
mass, scale radius, pericentre distance, and apocentre dis- 
tance with a sufficient grid of simulations, we use the com- 
putationally much cheaper rigid potential approximation for 
both, the main and the subcluster. 

Our simulations are run in the rest frame of the main 
cluster. The orbit of the subcluster through the main cluster 
is assumed to be that of a test particle falling through the 
main cluster, given an initial position and velocity. 

As the rest frame of the main cluster is no inertial frame, 
the ICM is subject to a pseudo-acceleration due to the at- 
traction of the main cluster core towards the approaching 



2.3 Model setup 

The Virgo cluster ICM is modelled by a spherically symmet- 
ric density and temperature profile, chosen to fit observa- 
tional data. Assuming hydrostatic equilibrium, we calculate 
the gravitational acceleration due to the underlying main 
cluster potential as a function of radius and thus the un- 
derlying cluster potential and mass profile, (see Sect 



2.4.1 



When the cluster is allowed to evolve for 2 Gyr in isolation, 
no sloshing or other modification of the ICM distribution is 
observed. 

Once we know the main cluster's potential, we calcu- 
late the orbit of a test mass moving through this potential. 
The subcluster will move along one of the orbits discussed 
in Sect. |2.4.3| Initially, the subcluster is placed at its orbit 
with a position and velocity such that it passes the pericen- 
tre after 1 Gyr. For the "slow", "medium" and "fast" orbits 
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Table 1. ICM parameters for the Virgo cluster model. 
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Figure 2. Overview of the Virgo cluster model, comparison with 
observations: Panels show electron density, temperature, pressure, 
and cumulative mass. For details see Sect. l2.4TT1 



(see Sect. 2.4.3) this corresponds to an initial separation of 
approximately 1 Mpc, 1.4 Mpc, and 1.6 Mpc, respectively. 
During the course of the simulation, the subcluster poten- 
tial is shifted along its orbit through the main cluster (see 
Fig.pl. The time normalisation is chosen such that pericen- 
tre passage happens at t = 0. 

We stop the simulations at t w 1.7 Gyr after core pas- 
sage when a good match with the observed CFs in Virgo is 
reached. This moment is well before the second core passage 
of the subcluster. 



2.4 Cluster models 

2.4-1 Main cluster - Virgo 

Figure [2] summarises the characteristics of our Virgo clus- 
ter model. The first panel shows observational data for the 
electron density. 



from the ACCEPT sample (Cavagnolo 
et al.|[2009 ) and the deprojected profile (Churazov, private 
comm.), along with a single- f3 model fit from Churazov et al. 



(2008), and a double- /3 model fit from Matsushita et al 



density: 


double (3 profile 


core radii 7-1,2/ kpc: 

core densities poi,02/(gcm — 3 ): 

/5i, 2 : 


1.20 23.68 

4.42 • 10- 25 1.38 • 10- 26 

0.42 0.52 


temperature: 


see Eqn. ill 



slopes 77ii 5 2 /(Kpc 1 ): 


-5.6 




0. 


Toi,o 2 /(10 7 K): 


1.78 




3.12 


break radius r*b r eak/kpc: 




27.7 




break range a<break/kpc: 




10 





( 2002 ) . In order to model the Virgo cluster as accurately as 



possible, we fit the deprojected data with a double-/^ model. 
The second panel shows the observed projected temperature 
profile from the ACCEPT sample, the deprojected data, and 
a fit from Churazov et al. ( 2008 ). We fit the deprojected data 
with 



T(r) 
Ti(r) 
T 2 (r) 

n 2 

C{r) 



(1) 



= 1 



Ti(r)C(r) + T 2 (r)[l - C(r)], where 
mi • r + Toi, 
m 2 • r + T 2 , 

(?71l - m 2 )r b reak + Til <4> Tl(r b reak) = T 2 (r b reak) 
1 



1 +exp (- r rbr ^ k ) ' 

\ a break / 



which connects two linear functions smoothly. At large radii, 
we use a constant temperature, contrary to the fit of |ChuraT] 
zov et al. (12008]) . The next panel shows the pressure profiles 



resulting from the ACCEPT data, the deprojected data, the 



Churazov et al. (2008) fits, and our fits to density and tem- 



perature. Assuming hydrostatic equilibrium, we can calcu- 
late the underlying potential and the associated cumulative 
DM mass, which is shown in the fourth panel, along with 
the cumulative mass for the ICM atmosphere. Table [I] sum- 
marises our ICM fit parameters. 

We study two different initial metal distributions, a 
steep and a natter one, which are described in Sect. |5J"| 



24.2 Subcluster 

The gravitational potential of our subclusters is described 



by a Hernquist halo (Hernquist 1990). We vary the mass 
and size of the disturber between 0.5 and 4 • 10 13 Mq and 
between 50 and 200 kpc, respectively (see Table [3j. At the 
upper end of our mass range, we are restricted by the con- 
dition that we want to model minor mergers only. A mass 
of 4 • 10 13 Mq equals the mass of the inner 300 kpc of Virgo. 
Hence, an interaction with higher mass subclusters would 
not be a minor merger any more and would leave obvious 
imprints in the Virgo cluster, which are not observed. The 
lower end of the mass range is determined such that the 
subcluster still leaves an imprint that is comparable to the 
observations. 

We abbreviate the combinations of subcluster mass, 
size and the pericentre and apocentre distance of its or- 
bit as MXaydminZdmaxQ, where "MX" shall indicate 
a subcluster mass of X • 10 13 Mq, aY a scale radius of 



6 E. Roediger et al. 



Y kpc, dminZ the pericentre distance of Zkpc, and dmaxQ 
the apocentre distance of QMpc. E.g. in our fiducial run, 
M2al00dminl00dmax3, a subcluster with 2 • 1O 13 M and 
scale radius 100 kpc passes the cluster centre with a mini- 
mal distance of 100 kpc, moving along an orbit where the 
maximal separation between both clusters was 3 Mpc. 



24.3 Orbits 

Slower subclusters on orbits with smaller pericentre dis- 
tances spend more time close to the cluster core and thus 
should cause a stronger impact. We test these expectations 
by modelling the passage of the subcluster along different 
orbits, whose characteristics are summarised in Fig. [3] We 
study orbits with pericentre distances, d m i n , of 20, 100, and 
400 kpc. We construct three groups of orbits, namely "slow" , 
"medium", and "fast, which arise from a maximum initial 
separation between subcluster and main cluster of 1, 3, and 
10 Mpc, respectively. 

The subcluster approaches the main cluster core from 
the — ^/-direction, passes left of the main core at t = 0, and 
recedes towards the +?/-direction. 



3 COMPARISON WITH OBSERVATIONS - 

INFERRING MERGER CHARACTERISTICS 

In our simulations, the gas sloshing proceeds as described in 
Sect. [TTT] and is qualitatively independent of subcluster and 
orbit characteristics. In order to focus the reader's attention 
to the aims of this work, we give a brief description of the 
sloshing dynamics in Appendix [B] 

With the aim of inferring the properties of the respon- 
sible subcluster and its orbit, we here present a detailed 
comparison between true and synthetic observations of the 
CFs in the Virgo cluster and apply our results also to other 
CF clusters. 

Qualitatively, the synthetic observations of all sim- 
ulation runs are similar. As the representative case, 
we show maps of X-ray brightness, brightness resid- 
ual^] and projected temperature for our fiducial run 
M2al00dminl00dmax3 in Figs. [4] [5] and [6] respectively. 

The preparation of our synthetic observations is de- 
scribed in Appendix [C] 



3.1 Basic evidence for sloshing scenario 

3.1.1 Brightness and temperature edges 

In the simulations, the gas sloshing produces observable arc- 
like X-ray brightness edges (Figs. [4] [5]) and analogous dis- 
continuities in temperature maps (Fig.pl. At these edges, 
the inner, brighter side is always the cooler one. The first 
edge appears towards the -\-x and +?/-direction of the clus- 
ter core about 0.2 Gyr after the pericentre passage (see also 
Fig. Bl). At the opposite side, a front at a smaller radius 
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becomes evident after 0.5 Gyr. With time, the edges move 
outwards. They can be seen along all LOSs. 

Thus, the basic characteristic of the sloshing CFs, 
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Figure 3. Different orbits of the subcluster: The top figure dis- 
plays the orbits in the xy-plane, one pericentre distance per panel. 
The line thickness codes the orbital velocity: slow, medium and 
fast orbits arise from an initial separation between subcluster and 
main cluster of 1, 3, and 10 Mpc, respectively. The circles along 
the orbit are spaced by 250 Myr. We initialise our subclusters in 
the "lower" (— y-direction) part of the orbital plane. The pericen- 
tre passage happens at t = 0. The bottom figure displays the 
evolution of the distance to the cluster centre, the total velocity, 
and the ram pressure as a function of time for each orbit (see 
legend, pericentre distances are given in kpc). 



1 residuals w.r.t. the azimuthal average 
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Figure 4. X-ray images: Synthetic images for fiducial run (high- 
resolution case, compare to Appendix [A| at t = 1.7 Gyr along 
three LOSs. Top left: Along the z-axis, which is perpendicular 
to the orbital plane. The subcluster orbit is marked by the black 
line. Bottom left: Along the y-axis, i.e., along the orbit. Top 
right: Along the x-axis. The brightness is given in arbitrary log- 
arithmic units (see colour scale). Bottom left: Observed X-ray 
image (same as in S10, but consistent colour scale), oriented such 
that NW is right and NE is up. 




Figure 6. Projected temperature maps corresponding to Fig.^] 



namely arc-like brightness edges accompanied by positive 
temperature jumps on opposite sides of the cluster core, are 
found in simulations and observations (Fig. [I]). 



X-ray brigrrtness/azimuthal average axis: z t=1700.35Myr 




Figure 5. Brightness residual maps corresponding to Fig.^] For 
the observed map, compare to Fig. ^ 



3.1.2 Spatially congruent cool structures and brightness 
excess structures 

In addition to the actual discontinuities, sloshing produces 
a characteristic cold pattern in projected temperature maps 
(Fig. [6] see also AM06). It takes the form of a single cool 
spiral winding outwards from the cluster core if the clus- 
ter is seen perpendicular to the orbital plane, and the form 
of cool arcs on opposite sides of the cluster centre if the 
LOS is parallel to the orbital plane. The arcs are always 
found at the axis belonging to the orbital plane, never at 
the perpendicular one. The outer edges of the cool spiral or 
cool arcs are the actual CFs. The brightness residual maps 
show exactly spatially congruent excess patterns (Fig. |5]). 
This correspondence of brightness excess and cool regions 
is found in all known sloshing CF clusters where both data 
sets exist. In the case of Virgo, the temperature and metal- 
licity maps (Simionescu et al. 2007) extend only out to 70 
kpc from the cluster centre, but clearly show the SE CF and 
its inner cold, metal-enhanced arc to be spatially congruent 
with the brightness excess feature in Fig. [I] 

The brightness excess spirals or arcs are the most easily 
detectable signatures of gas sloshing that are seen even with 
poor resolution, when the fronts themselves may be smeared 
out. The morphology of the excess structures may also be 
more easily detectable than the one of the edges themselves. 
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3.2 Morphology and orientation 

3.2.1 Our line-of-sight (LOS) is perpendicular to the 
subcluster 's orbital plane 

The morphology of the observed brightness residual maps 
presented in Fig. [I] depends to some degree on the method 
of its generation. A spiral morphology of both, the bright- 
ness excess region and the fronts themselves, is clearly seen 
in the residuals w.r.t. azimuthal symmetry, and is also sup- 
ported by the residuals w.r.t. the best-fit /3-model. Hence, 
from comparison to the simulations (Sect. 3.1.2 and Fig.|5|, 
we conclude that our LOS towards the Virgo cluster is ap- 
proximately perpendicular to the orbital plane of the per- 
turbing subcluster, and focus on this configuration for the 
further comparison. 

We note that we detect the CFs in projections along all 
grid axes, which was also shown by AM06. Hence, contrary 
to the suggestion of Birnboim et al. (2010), the detection of 
sloshing CFs in a given cluster does not require a favourable 
viewing angle. However, the signature is slightly weaker if 
the LOS is parallel to our x-axis, i.e. parallel to the orbital 
plane but perpendicular to the orbit. 



3.2.2 Orientation on the sky: NNW corresponds to our 
-\-x- direction 

In the Virgo cluster, the outermost CF is found NNW (la- 
belled "NW" for short) of the cluster core, the second front 
is on the opposite side (labelled "SE front"). The brightness 
excess spiral goes anti-clockwise from the centre outwards. 
We now rotate the observed brightness residual map such 
that the orientation of the observed brightness excess spiral 
matches the one in our synthetic maps (Fig. [5}. The best 
match, regarding the orientation of the excess spiral as well 
as the brightness excess distribution at the outer boundary 
of the observed field of view, is achieved by identifying the 
NW front with our H-x-direction or up to 45° clockwise from 
this. The same orientation is also favoured when comparing 
the radii of the NW and SE front between observation and 
simulation (see Sect. 3.3). For further comparison, we iden- 



tify the direction of the NW front with our x-axis. 

The current position of the perturbing subcluster is dis- 
cussed in Sect. [831 



3.3 The size of CF spiral constrains the age to 
- 1.5 Gyr 

In Fi g. [7| we measure the radii and widths of the CFs (see 
Sect. |C3| in all simulations into the SE, NW, and NE di- 
rection and plot their evolution as a function of time. For 
all subclusters and orbits (except M4al00dminl00dmax3), 
the CF radii depend mostly on the age (i.e. the time since 
pericentre passage) and very little on the characteristics of 
the subcluster and the orbit. There is a weak trend that 
configurations with less impact (small mass, large scale 
radius, large pericentre distance, fast passage) result in 
marginally smaller CF structures. Hence, we plot the re- 
sults for all simulations in Fig. [7] by grey lines and highlight 
only the fiducial run (M2al00dminl00dmax3) and the ex- 
ception M4al00dminl00dmax3 (discussed below). The radii 
differ between different subclusters and orbits only for the 



outermost CF in the NE (+y) direction, the direction the 
subcluster recedes to. However, going anti-clockwise from 
the +2/-axis, the sharp outer edge of the cool spiral turns 
into a flattening gradient in temperature (Figs. Bl and 16]), 
i.e. the true CF turns into a cool "fan". The radius of the 
outer boundary of the cool fan depends on the subcluster 
and the orbit characteristics, but also here we find a high 
degree of degeneracy. 

The only exceptional subcluster and orbit configuration 
is M4al00dminl00dmax3, where the rather massive subclus- 
ter is close to disrupting the cool core. Here, at a given time, 
the CF radii are significantly larger. 

Apparently, the sloshing itself is governed mainly by 
the underlying potential of the Virgo cluster. Thus, we can 
derive the age of the CFs in Virgo from the CF radii. In 
Fig. [7| we mark the radii of the observed Virgo CFs by hori- 
zontal black lines (inner and outer edge of each CF, see also 
Table [2ft. For an age of 1.7 Gyr, our simulations reproduce 
the size of the observed brightness excess spiral in NW and 
SE direction simultaneously. In the NE direction, most of 
the subclusters which do not reproduce the contrasts of the 
X-ray brightness and the temperature across the CFs cor- 
rectly (see Sect. 3.4 and Table |3| also produce the cool fan 
at an either too small or too large radius. The strongest im- 
pact subcluster M4al00dminl00dmax3 does not match the 
CF radii in the NW and SE simultaneously at any given 
time and any reasonable orientation. For the remaining con- 
figurations, alternative orientations of our synthetic maps 
on the sky fail to simultaneously match both observed CF 
radii, supporting our earlier constraint on the orientation 
(Sect. 



3.2.2). 



Given the simplicity of our model, the size of the 
observed CF structure is reproduced remarkably well at 
t — 1.7 Gyr. The direct comparison (Roediger & ZuHone, in 
prep.) between rigid potential simulations such as the ones 
used here and full hydro+Nbody simulations shows that the 
rigid potential simulations overestimate the age of CF by 
~ 200 Myr. Hence, we conclude that the CFs in the Virgo 
cluster are about 1.5 Gyr old. This more than doubles the 
estimate derived in S10. 



3.4 Contrasts across the edges constrain the 
subcluster mass/size/orbit 

With time, the CFs move outwards, and the contrasts of 
all quantities across them decrease. The age of the CFs is 
already fixed by their radius. Thus, we compare our simu- 
lation results at simulation time t = 1.7 Gyr (see Sect. 3.3) 
to the observations and use the contrasts to constrain the 
subcluster properties. In brief, subclusters with more mass, 
smaller scale radius, small pericentre and apocentre distance 
cause higher contrasts. 

From radial profiles like in Fig. [8] (see Appendix C2 for 
their calculation), we measure the projected temperature 
and X-ray brightness at the inner and outer edge of the 
NW and SE fronts, and derive the contrast of the X-ray 
brightness as the ratio of the values at the inner and outer 
edge (see Appendix C3). The results derived in this manner 
from the observation as well as the fiducial simulation are 
listed in Table 12 

Our automatic CF detection 



see Sect. C3) leads to 



smaller CF widths and, consequently, to lower X-ray con- 
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Figure 7. CF radii as a function of time. From temperature profiles towards the directions -\-x (NW), — x (SE), and -\-y (NE), we 
derive the radius (see Sect. |C3| of the outermost CF. For the NW and NE directions, we also show the radius of the second CF or 
cold fan, see labels in the middle and rhs panel. The thin black horizontal lines mark the radius of the observed CF structure. The thin 
vertical black line marks the time of 1.7 Gyr. We show the results of all simulations runs in light grey lines and highlight the fiducial run, 
M2al00dminl00dmax3 (see legend). The CF radii depend mostly on time and very little on subcluster and orbit characteristics, except 
for the cold fan in the NE. Only in the case where the most massive subcluster passes close to the Virgo core, M4al00dminl00dmax3 
(see legend), are the CF radii significantly larger than in all other cases. The error bars indicate the width of the fronts. 



Table 2. Quantitative characteristics of observed and simulated CFs (fiducial case). The observed values are derived from the profiles 
of S10. For the simulation, we give the values as derived by the automatic CF detection described in Appendix |C3| In brackets, we give 
the resulting X-ray brightness contrast and temperatures assuming the same radii as in the observation. 



quantity 



NW 
observation simulation 



SE 



observation 



simulation 



radius in kpc 80 - 100 83 - 95 (80 - 100) 32 - 50 40 - 48 (32 - 50) 

X-ray brightness jump 2.28-2.46 1.72(2.1) 2.11-2.18 1.51(2) 

Temperature inside CF in keV 2.29 - 2.34 2.24 (2.2) 2.4 - 2.43 2.43 (2.4) 

Temperature outside CF in keV 2.58 - 2.72 (< 3) 2.68 - 2.71 < 2.8 

Metallicity jump 1.15 - 2. 1.5 - 1.8 



trasts over the front, because part of the general decrease of 
brightness over the CF width is not included. If we assume 
for a CF width the same values as in the observation, the 
simulation nearly reproduces the observed jump in X-ray 
brightness. Also the temperatures just inside the CFs are 
reproduced well, whereas the over-prediction outside of the 



CFs is an artefact of the simulation method (Sect. 2.2). 

We derive the contrasts for all our simulations, com- 
pare them to the observed ones, and list the outcome of the 
comparison in Table [3] The runs in the bottom group in 
this table produce either too strong or too weak contrasts. 
Interestingly, these runs also predict a different radius of 
the CF/cool fan towards the NE (Fig. [7|), where too weak 
contrasts correspond to a too small radius and vice versa. 
Hence, we can reject these runs consistently from, both, the 
size and shape of the CF structure in maps and from the con- 
trasts across the CFs. For completeness, we compare their 
synthetic X-ray brightness and temperature profiles to the 
observed ones in Fig. |D1| 



3.4-1 Degeneracy of subcluster + orbit characteristics 

Not only our fiducial run, but seven of our runs (first two 
groups in Table pi reproduce the observed profiles nearly 
equally well. We plot projected radial profiles for X-ray 



brightness, temperature, and metallicity towards the NW 
and SE direction for all these runs in Fig. [8] 

3.4.1.1 Similarities Like the profiles, also the synthetic 
maps of all quantities differ remarkably little between these 
runs. We show brightness residual maps of all best-match 
cases in the appendix in Figs. |D2| and |D3| All of them re- 
produce the observed brightness residual map equally well. 

3.4.1.2 Differences Distinguishing observationally be- 
tween the best-match configurations is very difficult. 

A close inspection of the profiles reveals that in some 
cases (second group in Tablep] grey lines in Fig.pl there is a 
brightness deficit compared to the observation around radii 
of 20 or 30 kpc. This is due to a well- formed secondary CF 
towards the NW inside this radius, which is not detectable 
in the observed X-ray profiles. Also the remaining runs pre- 
dict this secondary CF, but exhibit only a small deviation 
from the observed X-ray profiles. Although the Virgo X-ray 
profile does not show this CF, both, the temperature and 
the metallicity exhibit a steep gradient in the NW profile at 
20 kpc. It is also evident in the metallicity map of Simionescu 



et al. 



(|2007|. This is exactly where the simulations predict 
a secondary CF in the NW (see also Fig. |7|. The obser- 
vational signature of this CF may have been disturbed to 
some degree by the nearby active galactic nucleus (AGN) 
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Table 3. List of simulation runs and their (dis) agreement with observations: Subcluster mass (column 1), subcluster scale radius 
(col. 2), pericentre distance of orbit (col. 3), apocentre distance of orbit (col. 4), result of two comparison criteria (col. 5, 6; exclamation 
marks indicate strong mismatch, brackets weak mismatch, "-" indicates no mismatch). The following comments (col. 7) apply to some 
simulations: (a) No simultaneous match of CF radii, close too disruption of core, (b) Perturbed X-ray profiles, (c) Too strong secondary 
CF in NW, brightness deficit inside 40 kpc in SE. 

The first and second group are our best-match simulations, configurations in the second group differ from the observations in details. 
The configurations in the third group do not fit the observation, remarkably the mismatch occurs consistently in more than one criterion. 
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in the Virgo centre. In the simulations, the differences re- 
garding the secondary NW CF are so subtle that we cannot 
disregard any of our seven best-match cases on this basis. 

The large-scale asymmetry described in more detail in 
Sect. |4.2| provides the most promising means to break the 
degeneracy: it is strongest for massive subclusters with large 
pericentre distances and weakest for the smallest pericentre 



distances (Figs D2 and D3). However, also here the differ- 
ences are subtle and may be hard to distinguish from intrin- 
sic asymmetries in Virgo. 



3.4-2 Comment on metallicity 

Given that the metallicity contrast across the CF depends 
strongly on the initial (unknown) metal distribution, it is 
not straightforward to use the observed contrast as a con- 
straint on sloshing dynamics. However, for our initially flat- 



ter metallicity profile (see Sect. 5.1 ), our simulated metallic- 



ity profiles match the observed ones inside 100 kpc (Fig. [8| 
well. Outside ~ 100 kpc, our simulations produce a signif- 
icantly lower metallicity compared to the observed values 
and, consequently, a too strong contrast across the NW 
front. This is due to an inaccurate initial metallicity at 
these large radii. Adding initially a metallicity floor outside 
~ 100 kpc would reduce the contrast across the NW front 
without changing the contrast of the SE front, which is at a 
smaller radius. As the sloshing increases the metallicity out- 
side ~ 100 kpc in the SE direction, this constant floor would 
be turned into a slightly rising metallicity profile, just as 
observed in this region. 



4 OBSERVABLE FEATURES OF EVOLVED 
CFS 

The basic observable properties of sloshing CFs as discussed 
in previous work and in the previous section are 

• spatially congruent cool, brightness excess features in 
form of a one-armed spiral or arcs on opposite sides of the 



core, depending on the LOS (Sect. 3.1.2) 

• brightness and temperature discontinuities, i.e. the ac- 
tual CFs, at the outer edges of the cool spirals or arcs 
(Sect." 



3.1.1) 



• steep positive gradients in temperature profiles and 
steep negative ones in X-ray brightness and metallicity at 



the radii of the CFs (Sect. 3.4) 



We here refine and extend this list and our understanding 
of the features. 



4.1 Alternating structure of cool+brightness 
excess and warm+brightness deficit 

Sloshing produces a characteristic staggered pattern: Re- 
gions of lower temperature and surface brightness excess 
alternate with brightness deficit regions of higher temper- 
ature, either in a spiral or staggered arcs pattern depending 
on the LOS. This pattern can be seen in maps (Figs. p\ 
[6l as well as in profiles (Fig. [8b. It even extends beyond 
the outermost CF. It is seen in all simulations and also in 
observa tions: in Virgo (Fig. [I), A496 ( |Dupke et al.||2007| ), 
A2142 flMarkevitch et al.pOOu), MS 1455.0+ 2232 and RX 
J172Q. 1+2638 flMazzotta fe Giacmtucci|2008 ). 

The cool regions are congruent with the brightness ex- 
cesses because the cool gas is associated with a higher gas 
density, which leads to a higher X-ray brightness and vice 
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Figure 8. Comparison of azimuthally averaged profiles of X-ray 
brightness, projected temperature, and projected metallicity 
for best runs. The thin vertical lines mark the observed CFs. 
The dashed (red) lines are for the cases that match also the 
overall X-ray brightness profiles best. The light grey lines collect 
the cases that match the profiles at the CFs but show small 
deviations in the inner X-ray brightness profiles. 
For the metallicity, we show the results for two initially different 
metal distributions: a steep one that follows the light profile of 
M87 (third row), and an initially flatter one (fourth row, see 
Sect.pTl). 
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Figure 9. Brightness residual maps for a large field of view, cor- 
responding to Fig. ^] For the observed map, compare to Fig. ^ 
The red circle in the top left panel indicates the size of the ob- 
served field of view. 



view. Except for the central spiral excess, there is no obvious 
trace of the interaction. The residual maps, however, (Fig.M) 
reveal that there is a clear asymmetry in the X-ray bright- 
ness. The central brightness excess spiral continues from the 
+?/-direction anti-clockwise, forming a brightness excess fan. 
The outer boundary of this fan is not a sharp edge. Relative 
to the orbit of the subcluster, the fan is located at the same 
side of the main cluster core as the subcluster's pericentre. 
Outside of the fan, there is a large brightness deficit region, 
whereas on the opposite side, there is a large brightness ex- 
cess region outside ~ 200 kpc. 

The large-scale asymmetry is strongest in the orbital 
plane, perpendicular to the cluster orbit. Also the metallicity 
is asymmetrical, where higher metallicity is associated with 
brightness excess. The asymmetry originates from the offset 
of the main cluster core w.r.t. the overall cluster towards the 
pericentre of the orbit. The same effect takes place in the 
full hydro+Nbody simulations (AM06, their Fig. 4). 

The detection of such a large-scale asymmetry could be 
used to distinguish the standard sloshing scenario from other 
scenarios, which could cause different large-scale signatures. 



versa. The same pattern can be deduced from the mock ob- 
servations of AM06. 



4.2 Asymmetry on large scales 

From our simulations, we also derive synthetic maps out to 
500 kpc from the cluster centre. At first glance, the large- 
scale X-ray images appear approximately symmetrical and 
undisturbed as soon as the subcluster has left the field of 



4-2.1 Large-scale asymmetry in Virgo 

The observed residual map of Virgo extends to about 170 
kpc. Inside this range, it agrees with the predicted asym- 
metry. It will be interesting to compare our prediction to 
larger-scale observations. The detection of such a large-scale 
feature would greatly help to constrain the orientation of 
our simulated maps on the sky. However, Virgo is a dynam- 
ical cluster, still in the state of assembly, so that also other 
processes will influence the symmetry or asymmetry at this 
length scale. 
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4-2.2 Large-scale asymmetry in A496, Perseus, 2 A 
0335+096 



The cluster A496 ( Tanaka et al.|2006 ) shows a clear bright- 
ness excess (and also cooler) region ~ 400 kpc towards N- 
NW, from which the authors derive that the subcluster 
passed the main cluster core from NW to SE. Our simula- 
tions indicate that the subcluster moved in W-E direction, 
passing the cluster core in the S. The orientation of the 
spiral formed by the CFs observed by Dupke et al. (2007) 
indicates for a subcluster motion from W to E. This excess 
region also matches the orientation of the central CFs de- 



scribed in Dupke et al. (2007), where the outermost CF is 
towards N. Additionally, Tanaka et al. (2006) report another 



CF about 200 kpc S-SE of the cluster centre, which is evi- 
dent as a brightness excess region in their residual map. Our 
simulations suggest that this is the cool brightness excess 
fan, which should be and is found opposite of the large-scale 
excess region. Thus, our simulations seem to reproduce also 
this cluster qualitatively. A496 will be an interesting target 
for future simulations aiming at a quantitative reproduction 
of the observed features. 

The residual and temperature maps for the Perseus 
cluster ( [Churazov et a l. 2003) qualitatively show exactly 
the same features: a cool brightness excess at large radius 
towards E, a large brightness deficit region at large radii 
towards W. The large brightness excess region at smaller 
radius W of the core could be the cool fan. North of the 
cluster core, there is a faint connection between the large- 
scale excess in the East and the fan in the W, which is also 
seen in our simulations. This argues for a subcluster orbit 
from N to S, passing W of the cluster core. We note that 
there is an inner brightness excess spiral in Perseus that 
does not seem to continue the outer structures, which could 
indicate a second sloshing event. 

The features in the large-scale images of the cluster 2A 
0335+096 ( Tanaka et al.|2006 ) show considerably more sub- 
structure, but are still consistent with a brightness excess 
fan towards the NW and a large-scale excess towards the 
SE. The fan could be the continuation of the central cool 



spiral seen by Sanders et al. (2009a). If this interpretation 



is correct, the subcluster should have moved from SW to- 
wards NE, passing the cluster core in the NW. Tanaka et al. 
(2006) discuss cooling or ram pressure stripping or turbu- 



lent viscous stripping as the origin of this cool brightness 
excess region. Our simulations suggest that they could be a 
by-product of the passage of a subcluster that was already 
stripped before pericentre passage. 



4.3 Substructure in profiles 

Radial profiles taken towards the CFs (Fig. p| show a con- 
sistent structure in both, simulations and observations: 

• When approaching the CF from the inside, the tem- 
perature decreases or forms a plateau. Also the metallicity 
shows a plateau. 

• The front itself is manifested as a steep increase in 
temperature and a steep decrease in metallicity and X- 
ray brightness. The apparent width of the CFs is due to 
azimuthal averaging over concentric ring segments whose 
boundaries do not coincide with the CF. The contrast of 



any quantity across a CF, defined as the ratio of the quan- 
tity at the inner and outer edge of the CF, always includes 
the intrinsic variation over this radial range. 

• Outside the CF, the temperature again decreases, 
i.e. there is a warmer rim enclosing each CF. 

The same structure is found in other CF clusters, as far 
as data quality allows its detection. The decrease of tem- 
perature with radius inside the CFs may be a speciality of 
the compactness of the Virgo cluster, but plateaus in tem- 
perature inside the CF are found for the N front in A496 



(Dupke et al. 2007), the SW front 2A 0335+096 (Sanders 



et al. 2009a), in A1795 (Markevitch et al. 2001), A2142 



flMarkevitch et al.||2000|), A2204 (|Sanders et al.||2005[ ), MS 



et al. 



et al 



1455.0+2232 (Ma zzotta fc Giacintu cci 2008), NE f ront in 
NGC5098 JRandall et al.||2009| ), probably in A120 1 flOwers 



2009). A metallicity plateau is seen in A2204 \ Sanders 
2009b). Indications of a warm rim outside the CF are 



seen in A1795, A2204 and A2142. However, most published 
data does not resolve the temperature profiles well enough 
to clearly detect these features. Occasionally, only the tem- 
perature jump across the CF is measured. 

For future observations we suggest to derive the full 
profiles from the cluster centre beyond the CFs. This will 
enable a better comparison to simulations, and tighter con- 
straints on the identification of the responsible subcluster. 
Also an alignment of the spectral extraction rings with the 
CFs could be helpful in measuring accurate contrasts. 



5 METAL TRANSPORT 

While sloshing redistributes the central gas, it must also 
redistribute the heavy elements throughout the central re- 
gion. In addition to producing metallicity discontinuities at 
the CFs, which are observed in Virgo (Fig.pl, sloshing may 
also be able to broaden the metal distribution throughout 
the cluster centre. We study this process by tracing the evo- 
lution of an initially centrally peaked metallicity distribution 
in our fiducial run. We find the evolution of the metal distri- 
bution to be tightly coupled to the overall dynamics. Sim- 
ulations that produce similar temperature structures, also 
produce similar metallicity structures. 



5.1 Initial metal distribution 

We study an initially steep and a natter radial metal density 



profile. In the steep case, we follow |Rebusco et al. (2006), 
who assume that originally the central metal distribution 
in a cluster follows the light profile of the central galaxy 
(here M87, |Kormendy et al.|2009| ), but is broadened by dif- 
fusion caused by turbulent gas motions. It is also debat- 
able whether the metal distribution about 2 Gyr ago should 
follow the current light profile of M87. If the metal peak 
in the Virgo centre was produced by stellar mass loss and 
SNIa, then the O/Fe profile should not be fiat as observed 
(S10). More likely, the metal distribution originates from 
the assembly of M87 and has already undergone some re- 
distribution. Hence, we also consider a flatter metal profile 
which approximates the observed metal density inside the 
NW CF. In both cases, we fit the initial metal density by 
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Table 4. Parameters for deprojected Sersic profiles describing 
the steep (fitted to M87 light) and the flatter metal distribution. 

steep (oc M87 light) flatter 



scale radius 
Sersic index 



54.56 kpc 
11.84 



120 kpc 
2 



log(Emission-weighted metallicity / solar) axis: z 



log(Emission-weighted metallicity / solar) axis: z t=1700.06Myr 
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Figure 10. Synthetic projected metallicity map (in solar units) 
for initial state (left) and 1.7 Gyr after subcluster pericentre pas- 
sage (right), for the case of the flatter initial metal distribution. 



a deprojected Sersic profile ( Sersic|1963 Prugniel & Simien 
1997) with the parameters listed in Table 14] 

The evolution of 3D metallicity profiles across the CFs is 
shown in Fig. |B2| Metallicity discontinuities appear exactly 
at the CFs, but the contrast across them is clearly higher 
for the initially steep profile. 



5.2 Comparison to observations 



Projected metallicity maps (Fig. 10) display a spiral struc- 



ture coincident with the one seen in the temperature and 
brightness residual maps (Figs. pip]). 

The two bottom panels in both subfigures of Fig.[8]com- 
pare the simulated to the observed metallicity profiles. The 
steep gradients at all three CFs are clearly reproduced, so 
are the plateaus inside the CFs and the NW-SE asymme- 
try (see Sect. El. In case of the initially steep profile, the 
metallicity contrasts across the fronts are higher than ob- 
served, and the metallicity in the central ~ 50 kpc remains 
too steep compared to the observations. This shows that the 
redistribution by sloshing alone is not sufficient to transform 
a steep metal peak originating from stellar mass loss of the 
central galaxy into the observed flattened one. Additional 
processes are required (see e.g 
and discussion in Sect 



5.3.2). 



Rebusco et al. 2005 2006 



For the initially flatter case, we achieve a good agree- 



ment between observations and simulations (see Sect. 3.4). 
Therefore, we will concentrate on this case for the further 
analysis. 
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Figure 11. Effect of sloshing on the metal distribution. Top: 
cumulative metal mass for initial time and t = 0.8, 1.7 Gyr. Sec- 
ond panel: Relative change of cumulative metal mass: profiles 
from previous panel divided by initial distribution. Third: spher- 
ically averaged metal density (thick lines) for the initial state and 
t = 0.8, 1.7 Gyr. The thin lines show the metal density along the 
+a>axis of the simulation grid. Fourth: Relative change of metal 
density profiles: profiles from the previous panel divided by the 
initial profile. Fifth: diffusion constant derived from the profiles 
in the third panel. 
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5.3 Broadening of the global metal distribution 

5.3.1 General 

We derive the cumulative metal mass and the spherically 
averaged metal density for different timesteps and compare 



them to the initial state (Fig. 11). Sloshing does not only 
influence the metal distribution towards the CFs, but broad- 
ens the overall metal distribution. Within 1.7 Gyr after 
pericentre passage, the central metal density is depleted by 
about 25%, while the metal density between 75 and 95 kpc 
increases by about 10 to 20%. This is somewhat less than 
the 20 to 30% derived in S10, but there only the difference 
between the NW and the SE sector was compared, while we 
considered the global broadening. Sloshing is, however, not 
a symmetric process, and the metal redistribution is much 
stronger if considered along a line directly towards a CF. 
Along the x-axis of our simulation box, the metal density is 
enhanced by a factor of 1.5 to 2 within the 20 kpc inside of 



the CF (third and 4th panel of Fig. 11 ) 



5.3.2 Diffusion description 

Assuming the heavy elements in a cluster core originate 
mainly from stellar mass loss of the central galaxy and are re- 
distributed by turbulence, Rebu sco et al.| (2005, 2006) derive 
an effective diffusion constant of the order of 10 cm 2 /s for 



several clusters, and 0.8-10 cm /s for Virgo. Graham et al 



(2006) find a value of 10 cm /s for the central region of 



the Centaurus cluster, but considerably smaller values out- 
side 25 kpc. Roediger et al.| ( |2007[ ) have studied the metal 
mixing induced by AGN- inflated buoyantly rising bubbles, 
which drag metal-rich gas from the centre to larger radii. 
They derived effective diffusion coefficients that depended 
strongly on cluster radius, ranging from 10 29 cm 2 /s at about 
10 kpc from the cluster centre to about 10 26 cm 2 /s at 50 kpc 
radius and 10 22 cm 2 /s at 100 kpc radius. 

In order to quantify the effectiveness of sloshing in metal 
redistribution, we calculate the effective diffusion coefficients 
from our simulations by comparing the spherically averaged 
density profiles at t — 0.8, 1.7 Gyr to the initial state (see 



Roediger et al. (2007) for details). The result is shown by 
the thick lines in the bottom panel of Fig. |11| It makes clear 
that diffusion is not a good approximation here, because 
we derive negative, hence unphysical, diffusion constants in 
some parts. This occurs e.g. inside the major CF, i.e. at those 
radii where sloshing deposits the metals it has removed from 
elsewhere. True diffusion would act differently: metals moved 
from the centre to some larger radius would not remain there 
but continue to diffuse outwards until a flat distribution is 
achieved. In contrast, sloshing removes metals from different 
radii and "piles them up" inside the CFs. 

In the inner part we derive a positive diffusion constant 
of the order of 10 26 cm 2 /s between 10 and 40 kpc, which 
is less than what was derived for AGN activity. This is not 
surprising, because the AGN sits in the very centre of the 
cluster and creates the strongest impact there. In contrast, 
the very central gas will be affected least by sloshing be- 
cause it is held more tightly by the deep central gravita- 
tional potential. Thus, AGN activity and gas sloshing are 
complementary processes, the former being more effective 
in the very cluster centre, the latter at larger radii. In a re- 
alistic case, a combination of both could be very effective. 



Furthermore, the metal mixing by AGN induced buoyant 
bubbles is, on average, a continuous process, because a se- 
ries of bubbles is inflated by intermittent AGN activity. The 
sloshing studied here happens in the aftermath of a single 
event. Hence, the diffusion efficiency averaged over the whole 
event is rather low. If the effective diffusion constant is av- 
eraged over a much shorter time, e.g. 50 to 200 Myr, it is 
at least a factor of 10 higher. Also repeated minor merger 
events may enhance the efficiency of metal redistribution. 



6 SUBCLUSTERS WITH GASEOUS 
ATMOSPHERES 

Previous work (AM06) stresses that the subcluster needs 
to be gas-free to produce sloshing CFs but leave the over- 
all cluster appear relaxed. Having constrained the mass and 
size of the responsible subcluster, we want to test if a suf- 
ficient subcluster can be completely ram pressure stripped 
significantly prior to pericentre passage. 



6.1 Initial setup 

We fill the fiducial subcluster M2al00 either with an ex- 
tended or with a compact atmosphere. For the subcluster 
atmosphere we assume a constant temperature with a cen- 
tral drop described by 



T(r) 



1 + (r/rro) 3 
D T + (r/rro) 3 ' 



(2) 



where To is the temperature of the subcluster gas at large 
radii, Tq/Dt the central temperature, and rro the scale ra- 
dius of the cooler core. Demanding hydrostatic equilibrium, 
we can calculate the density and pressure profile for the sub- 
cluster gas. For the simplest case of a constant temperature, 
T(r) = To, these are 



P(r) 
Po 



Po 



exp 



mpGM suh 
kT 



1 



r + a sub 



1 



(3) 



where m p is the mean particle mass in the subcluster atmo- 
sphere, and M su b and a su b are the mass and scale radius of 
the Hernquist potential describing the subcluster DM distri- 
bution. The slope in density and pressure is very sensitive to 
the temperature of the subcluster gas. The subcluster atmo- 
sphere truncates where its pressure drops below the ambient 
pressure of the main cluster's ICM. This truncation condi- 
tion determines the size and total gas mass of the subcluster. 
We adopt an ambient ICM pressure of 5-10 -13 erg cm -3 , typ- 
ical for the Virgo cluster at 1 Mpc (Fig. |2J). Then we chose 
the temperature and density of the subcluster gas such that 
the subcluster gas pressure drops below the ambient ICM 
pressure at a radius such that, in both cases, the subcluster 
contains about the same amount of gas in this pressure limit. 
Table [5] summarises the parameters for both atmospheres. 

6.2 Result 

On our fiducial orbit, the compact atmosphere is not com- 
pletely ram pressure stripped. On the contrary, its dense gas 
interacts strongly with the Virgo core gas and destroys the 
cool core completely. 

The extended atmosphere is ram-pressure-stripped just 
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Table 5. Parameters for the gaseous subcluster atmospheres (see 
Sect. |6.l] for details). 
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Figure 12. Fiducial orbit, fiducial subcluster M2al00, filled with 
extended gaseous atmosphere: The residual map and projected 
temperature map at t = 1.7 Gyr are similar to the case without 
gas (Fig. [5J, although more irregular. 



at the pericentre passage. However, its less dense gas can- 
not penetrate the central 100 kpc of the Virgo core, where 
sloshing can go on in a similar manner as in the gas-free 
cases. The resulting CF structures are also similar (Fig. 12), 



although more irregular. The CFs on the NW and SE sides 
are found at the same positions as in the gas-free case. The 
X-ray brightness profiles are, however, more distorted than 
in the gas-free cases and in the observation, reflecting the 
more irregular CF structure. The ability of the gaseous at- 
mosphere to disturb the sloshing is sensitive to the balance 
between the gas density and the ram pressure the subcluster 
experiences. For example, a subcluster with a half as dense 
atmosphere is already ram pressure stripped 400 kpc before 
the pericentre passage, and the central sloshing proceeds as 
in the initially gas- free case. 

The regularity of the observed brightness excess struc- 
ture suggests that the subcluster responsible for the CFs in 
the Virgo cluster has lost its gas well before its pericentre 
passage. Alternatively, the subcluster may have passed the 
cluster core at a large pericentre distance, such that the sub- 
cluster atmosphere does not penetrate the main cluster core, 
where it would erase the sloshing signature. AM06 present 
such a case with a mass ratio of 5:1 and a pericentre dis- 
tance of about 400 kpc. The CF and cool spiral observed in 



ratio 2:1) requires a rather large pericentre distance to pre- 
vent the disruption of the cool cores. In both cases, either 
the early loss of the gaseous atmosphere or the passage at 
a large distance, traces of stripped subcluster gas might be 
found at larger distances from the cluster centre. 



7 SLOSHING TRIGGERED BY A GALAXY'S 
BOW SHOCK 

Any process that offsets the ICM in a "bulk- like" manner 
will cause subsequent sloshing and formation of CFs. An al- 
ternative to displacing the ICM by gravitational interaction 
with a subcluster is the passage of a shock, proposed by |Chu-| 
razov et al. (2003). These authors presented a simulation of 



an idealised plane shock passing a hydrostatic cluster. 

We want to advance the idea of a shock as the cause of 
sloshing CFs and study the impact of a bow shock of a mas- 
sive galaxy. Such a scenario is also attractive when regarding 
merger rates: passages of a 2 • 10 12 Mq galaxy are more fre- 
quent than passages of > 10 13 Mq sub clusters. Thus, we let 
a galaxy pass the Virgo core at a distance of 400 kpc (or- 
bit "dmin400 fast" in Fig [3b. The large pericentre distance 
ensures that its ram-pressure stripped tail does not appear 
in images of the central region. A high velocity is needed 
to produce a sufficiently strong shock. The galaxy potential 
is described by a Hernquist halo with 2 • 10 12 Mq and scale 
radius 30 kpc. Initially, the galaxy contains a gaseous halo 
of 3 • 10 6 K and central density 10 _24 g cm -3 . The exact dis- 
tribution of the gas in the galaxy does not matter as long 
as the galaxy can retain some gas and hence a bow shock 



up to pericentre passage (see Fig. 13). The galaxy may also 
be a spiral galaxy. The bow shock passes the Virgo core 
about 200 Myr after the closest approach of the galaxy. The 
shock normal points towards the diagonal between the +x- 
and H-y-directions, which is also the main direction of the 
triggered sloshing. Fig. [13] also displays a brightness resid- 
ual map at 1.2 Gyr after shock passage, i.e. 1.4 Gyr af- 
ter pericentre passage of the galaxy. In the inner 120 kpc, 
the synthetic maps resemble the ones of the classic scenario 
(Fig. [5I except that they appear to be rotated by about 
45° ant i- clockwise. In order to match the observed map, we 
identify the direction of the NW CF with the diagonal be- 
tween the +x- and +?/-directions. Figure [14] compares the 
observed and synthetic profiles of X-ray brightness and pro- 
jected temperature at this moment. The radii of all CFs are 
matched simultaneously, and the contrasts across the NW 
front agree well. The simulated ones across the SE front are 
too strong. Also the X-ray profile on the NW side slightly 
deviates from the observed one. However, these details will 
also depend on the details of the interaction with the pass- 
ing galaxy. Thus, we conclude that also the bow shock of a 
galaxy passing the Virgo core as explained above could have 
caused the observed CF structures in the Virgo cluster. 

The best way to distinguish between this and the classic 
scenario is provided by the larger-scale structure: Both, the 
cool brightness excess fan between about 100 and 170 kpc 
SE of the cluster core, and the pronounced brightness excess 



Abell 1644 (Johnson et al. 2010) seems to be such a case. 
However, the fact that Abell 1644 is a major merger (mass 



outside 200 kpc NW of the Virgo core (see Sect. 4.2) would 
be absent in the galaxy bow-shock scenario. In this case, a 
large-scale view may reveal traces of the gas stripped from 
the galaxy. Where exactly this may be found would also 
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Figure 13. Sloshing triggered by a galaxy's bow shock. Top: 
Large-scale X-ray image at the moment when the bow shock 
crosses the cluster centre. The ram pressure stripped galaxy can 
be seen at the lhs. Bottom: Brightness residual map out to 250 
kpc at 1.4 Gyr after the galaxy's pericentre passage. 
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Figure 14. Projected and azimuthally averaged profiles across 
the CFs after bow shock passage: X-ray brightness, projected 
temperature (compare to Fig. [8|. We match the orientation of 
our simulated maps on the sky such that the direction of the 
observed NW CF is towards the diagonal between the -\-x and 
+y-axes. 



depend on the original gas distribution inside the galaxy 
and thus the galaxy's ram pressure stripping history. 

A possible candidate galaxy is discussed in Sect. |8.2[ 



8 DISCUSSION 

8.1 Reliability of our model 

The strongest concern about our model is the restriction to 
rigid potentials. However, in a direct comparison to full hy- 
dro+Nbody simulations, Roediger&ZuHone (in prep.) show 
that the rigid potential approximation reproduces the ori- 
entation of the sloshing signature correctly. From the size of 
the CF structure, the rigid potential approximation overes- 
timates the age of the minor merger slightly, by about 200 
Myr. The temperature inside the CFs as well as the den- 
sity inside and outside of the CFs is reproduced well. In 
our comparison to the Virgo data, we have used only these 
well-reproduced characteristics. 



8.2 Identifying the responsible subcluster 

According to the merger geometry derived in Sect. |3.2| the 
subcluster that perturbed the Virgo core has moved from 
~W-SW to E-NE, passing the Virgo core in the South about 
1.5 Gyr ago. Consequently, the current position of the sub- 
cluster (see Fig. [3J) corresponds to about 1.5-2.5 Mpc E-NE 
of the Virgo centre. The uncertainty in distance is due to 
the ill-constrained velocity of the subcluster. Our test parti- 
cle orbits do not take into account dynamical friction, which 
will cause the subcluster to be at a somewhat closer distance 
than derived from Fig. [3] Thus, our simulations suggest the 
perturbing subcluster about 1-2 Mpc eastward of the Virgo 
centre. The subcluster should still be moving away from the 
Virgo centre. 

The region in question is just outside the Virgo cluster 
catalogue ( |Binggeli et al.||1985} |1987[ |1993| ) . An inspection 
of sky maps reveals the following candidates: M60 (radial 
velocitMj 1117kms _1 , Virgo mean is 957kms _1 ), which is 
situated about 1 Mpc E-SE of the Virgo centre; further north 
from here are NGC 4654 (radial velocity 1046 km s" 1 ), NGC 
4639 (radial velocity 1018 km s" 1 ), NGC 4659 (radial veloc- 
ity 480 km s" 1 ), and NGC 4689 (radial velocity 1616 km s" 1 ) 
still further North. All of these except NGC 4659 have a ra- 
dial velocity close to the Virgo cluster mean, suggesting that 
their orbits are close to the plane of the sky. The massive 
elliptical galaxy M60 contains about 10 13 Mq within 76 kpc 
(Shen & Geb hardt||2010 ), which is close to the mass we re- 
quire. However, it still contains a gaseous atmosphere of its 
own. It might have passed the Virgo core at a large enough 
distance to prevent its ram pressure stripped gas to show up 
in the central part. In this case, its bow shock could have 
assisted the gravitational triggering of the sloshing. Alterna- 
tively, it may have acquired a new gaseous halo within the 
last 1.5 Gyr. 

If the CF signatures in the Virgo cluster were caused 
by the bow shock of a galaxy (see Sect. [7|), the disturbing 
galaxy should now be found about 1.5-2 Mpc N-NE of the 



2 Heliocentric radial velocities are taken from the NASA/IPAC 
Extragalactic Database (NED). 
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Virgo core. A possible candidate is the peculiar SO galaxy 
M85, about 2 Mpc N of Virgo. Its heliocentric radial velocity 
is 730 km s -1 , which is again close to the Virgo mean. This 
galaxy even shows some recent star formation in its centre 
that may have been triggered by the interaction with the 
ICM, but may also be due to interaction with its neighbour 



NGC 4394 (see Sansom et al. 2006 and references therein). 



8.3 Stability of CFs 

CFs are expected to be associated with shear motions, which 
should lead to the formation of Kelvin-Helmholtz instabili- 
ties (KHIs) and possibly to the destruction of CFs. We in- 
deed observe this instability in the later stages (after 1 Gyr) 
in our simulations, although it is not omnipresent. More- 
over, we find that the CFs are never destroyed completely, 
but are reformed continuously. Observing the KHI might be 
intrinsically difficult: in Appendix[A]we show that in projec- 
tion the CFs appear less disturbed in the case of the highest 
resolution, although here the KHI is clearly present. With 
such high resolution (CFs resolved with 0.5 kpc), the KHI 
takes place at small scales, which are averaged out by pro- 
jection. Thus, CFs may appear sharp in projection although 
the KHI is taking place. 



8.4 Impact of sloshing on the radio lobes 

The central AGN of the Virgo cluster has induced a rich 
structure inside the inner ~ 30 kpc of the ICM, including 
radio lobes and cavities. The inner radio lobes show a pe- 
culiar morphology: while the eastern arm extends straight 
out to about 20 kpc, the western arm is bent sharply south- 
wards at about 15 kpc from the cluster centre. S10 discussed 
the possibility that the western arm is bent due to the same 
sloshing motions that cause the CFs. Furthermore, the outer 
radio lobes are rotated about 50° clockwise w.r.t. the inner 
lobes. In Fig.[l5]we compare the central velocity fields from 
three of our best-match simulations to the observed radio 
lobes. In our images, N-NW is along the +x-direction, and 
we have rotated the observed radio contours accordingly. In 
the central 15 kpc, the velocity field is weak and unlikely to 
modify the AGN outflow. Outside 15 kpc, the velocity field 
exhibits a clear clockwise rotational pattern, but the details 
depend on the subcluster and the orbit characteristics. At 
the western side, the flow moves into the same direction as 
the bending of the jet arm, but the reason of the sharp bend 
is not obvious. However, the rotational flow pattern could 
have rotated the outer, and older, radio lobes to their cur- 
rent orientation within ~ 100 Myr, which is comparable to 



their ages derived by Owen et al. (2000) and Werner et al 



(2010). The latter authors also argue that both, the inner 



and outer radio lobes extend considerably along our LOS. 
Thus, the radio lobes might be affected be the sloshing only 
mildly. The details of the interplay of the AGN activity and 
gas sloshing have to be studied with a set of dedicated sim- 
ulations. 



9 SUMMARY 

We perform a set of hydrodynamical simulations of minor- 
merger induced gas sloshing and the subsequent formation of 



cold fronts. For the first time, these simulations aim at repro- 
ducing the characteristics of the cold fronts and associated 
structures not only qualitatively, but also quantitatively. In 
this paper, our target is the Virgo cluster. 

By comparing synthetic and real observations in great 
detail, we show that the sloshing scenario can reproduce 
the known characteristics: a pair of cold fronts, their radius, 
the contrast in X-ray brightness, projected temperature, and 
metallicity across the fronts, and a spiral-shaped brightness 
excess. The comparison suggests the presence of a third cold 
front 20 kpc N-NW of the cluster centre, whose detection 
was unsure so far. 

Furthermore, we identify additional, so far unreported, 
features typical for sloshing cold fronts, (i) X-ray bright- 
ness excess regions alternate with brightness deficit regions. 
Exactly spatially congruent structures are found in tem- 
perature and metallicity maps, where brightness excesses 
correspond to cool, metallicity enhanced regions and vice 
versa. This alternating behaviour is also evident when over- 
plotting profiles across opposite cold fronts, (ii) Inside each 
cold front, the temperature profile is constant or radially 
decreasing, and it is accompanied by a plateau in metallic- 
ity. (iii) On the outside, the cold fronts are bordered by a 
warm rim. (iv) On a larger scale of a few 100 kpc, there 
is a typical large-scale brightness asymmetry, which will be 
helpful in constraining the orbit disturbing subcluster. We 
can trace these new features not only in Virgo, but also in 
other clusters exhibiting sloshing cold fronts. 

Our constraints on the disturbing subcluster are the fol- 
lowing: From the size of the observed cold front structure, 
we constrain the age of the cold front structure in Virgo to 
about 1.5 Gyr. The presence of the spiral-shaped brightness 
excess suggests our line-of-sight to be nearly perpendicular 
to the orbital plane. The contrasts of X-ray brightness and 
temperature across the fronts, exclude subcluster masses be- 
low 10 13 Mq and above ~ 4 • 10 13 Mq. The subcluster mass, 
its size and the pericentre distance of the orbit are degen- 
erate. The limits of the possible range are (2 • 10 13 Mq, 100 
kpc pericentre distance), and (4 • 10 13 Mq, 400 kpc pericen- 
tre distance). The disturbing subcluster must currently be 
located east of the Virgo centre. Its orbital velocity is no 
well-determined parameter, thus its current distance can be 
constrained between 1 and 2 Mpc only. We suggest M60, 
currently about 1 Mpc east of the Virgo centre as the most 
likely culprit. 

Additionally, we quantify the metal redistribution 
caused by sloshing in Virgo. We demonstrate that sloshing 
alone cannot transform a peaked metal distribution propor- 
tional to the stellar mass density in M87 to the observed 
flattened one. We also demonstrate that diffusion is inept to 
describe this particular redistribution process. 

To complete our analysis, firstly we show that the dis- 
turbing subcluster could indeed be completely ram pres- 
sure stripped before pericentre passage as required by the 
standard scenario, and thus leaves no traces but the ob- 
served ones in Virgo. However, a compact gaseous atmo- 
sphere might survive beyond the pericentre passage and 
leave additional features in the Virgo ICM. Secondly, we 
demonstrate that the bow shock of a fast galaxy passing the 
Virgo cluster at ~ 400 kpc distance also causes sloshing and 
leads to very similar cold front structures. The disturbing 
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Figure 15. Zoom-in of the velocity field (white arrows) and density field (colour-coded in gem -3 ) at t = 1.7 Gyr in the orbital plane of 
the subcluster. We show of our best-match cases: Mla50dminl00dmax3 (left), M2al00dminl00dmax3 (middle), M4al00dmin400dmax3 
(right). In the middle panel, we overlay the observed radio lobes (|Owen et al.|2000). 



galaxy would be located about 2 Mpc north of the Virgo 
centre, M85 being a possible candidate. 
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APPENDIX A: RESOLUTION TEST 

The formation mechanism of CFs is accompanied by shear 
motions which are prone to KH instability. In this instability, 
the small-scale modes grow fastest, hence, in simulations, 
their evolution can be resolution-dependent. Furthermore, 
the Virgo cluster potential is very steep in the inner few kpc, 
which is not resolved sufficiently in our standard-resolution 
simulations. This results in a warm ICM core in the inner 5 
kpc after the passage of the subcluster, which is absent in 
the high-resolution cases. However, this very central region 
does not influence our results. 

We test our fiducial run in 3 resolutions. Figure |Al] dis- 
plays temperature slices in the orbital plane for our standard 
resolution, a higher resolution case where the resolution is at 
least a factor 2 better at each position, and a very high res- 
olution which is again a factor of 2 better everywhere. The 
overall structure, regarding size and contrast across CFs, 
is independent of resolution. The fronts are discontinuities 
within the resolution for all resolutions. Even with higher 
resolution, clear KH rolls appear only after 1 Gyr after core 
passage. At no time the CFs are completely disrupted. With 
better resolution, KHI takes place at smaller scales. The 
sloshing motion continuously reforms the CFs. 

In projection, the KHIs are less pronounced (see X- 
ray images, residual images, projected temperature maps 



in Fig. A2), for the very high resolution, they are hardly 



detectable. The large-scale signature reported in Sect. |4.2| 
does not depend on resolution. Neither do the azimuthally 
averaged profiles of X-ray brightness, projected temperature 



and metallicity (Fig. A3) 



APPENDIX B: EVOLUTION AND DYNAMICS 
OF SLOSHING IN VIRGO 



We demonstrate the evolution of the gas sloshing in Fig. |Bl| 
which displays the ICM temperature (colour coded) and ve- 
locity field (white arrows) in the orbital plane at different 
timesteps. The sloshing proceeds as described in Sect. |1.1| 
The only effect not mentioned there is the adiabatic heating 
of the ICM near the subcluster due to infall of ICM into the 



subcluster potential (see top panel in Fig. Bl). 

Qualitatively, the evolution is independent of the sub- 
cluster and orbit characteristics, which mainly influence the 
intensity of the contrast of density and temperature across 
the CFs (see Sect. [3}. 

In Fig. |B2| we take a closer look at the structure of the 
cool spiral and CFs by plotting temperature, density and 
pressure profiles along the x-axis. The CFs are indeed dis- 
continuities in temperature and density within our resolu- 
tion, but the pressure is continuous. However, the pressure 
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Figure Al. Impact of resolution in fiducial run: Colour-coded temperature slices in the orbital plane. White boxes demonstrate the 
resolution, one box corresponds to one "block" of 16 3 cells. Three different resolutions are shown: Left: standard resolution, innermost 
blocks have a resolution of 1 kpc. Middle: High resolution, innermost blocks have a resolution of 0.5 kpc. Right: Very high resolution, 
same block distribution as middle panel, but each block contains 32 3 grid cells, and the innermost blocks have a resolution of 0.25 kpc. 
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Figure A2. Synthetic brightness residual maps corresponding to Fig |Al| 



gradient shows a discontinuity just inside of each CF. In- 
side of each CF, there is a plateau in density accompanied 
by an increasing temperature towards smaller radii because 
the cool gas from the centre has been moved outwards. Go- 
ing further inwards, in many profiles the next CF can be 
detected. The metal profiles will be discussed in detail in 
Sect.H 

The typical velocities induced by sloshing in our fidu- 
cial model are 300kms _1 , which is about half the circular 
velocity inside the inner 30 kpc of the Virgo cluster. This 
supports the idea of Keshet et al. (2009), that the circular 



component of the sloshing motions contributes a significant 
radial acceleration working on the sloshing gas in the form 
of centripetal acceleration. 



APPENDIX C: PREPARATION OF 
SYNTHETIC OBSERVATIONS 

In order to highlight the simulation results, e.g., regarding 
the width of the CFs, we do not convolve our maps with an 



instrument kernel, nor do we reduce the resolution of our 
maps to match current observations. 



CI Maps 

• Synthetic X-ray images are derived by projecting n 2 VT 
(where n is the particle density and T the gas temperature) 
along the LOS. This results in unnormalised brightness im- 



• We calculate brightness residual maps by dividing each 
X-ray image by its azimuthal average. We favour this version 
over residuals w.r.t. a fitted /3-model because we want to 
highlight the deviation from the initial symmetry. 

• Reducing the 3D temperature structure of a galaxy 
cluster to a 2D map requires (weighted) averaging along the 
LOS for each pixel of the temperature map. A standard way 
is to calculate the emission-weighted projected temperature: 



fWTdz 
JWdz 



with W 



(Cl) 



Mazzotta et al. (2004) show that using the weights 
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Figure A3. Projected, azimuthally averaged profiles of X- 
ray brightness and projected temperature (like in Fig. [51, at 
t = 0.5, 1.2, 1.7 Gyr (see legend), for two resolutions (thick lines 
for high resolution, thin lines for standard resolution, compare 
Fig |Al]) . 
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: n/T 



3/4 



(C2) 



results in estimates for projected temperatures which are 
closer to the value derived by spectroscopical temperature 
fitting. Therefore we use this method. Emission- weighting 
results into very similar maps, except that the coolest re- 
gions appear slightly less cool, similar to what is found by 
Mazzottaet al.N2004ft. 



• In analogy to temperature maps, we derive metallic- 
ity maps by calculating the emission- weighted mean of the 
metallicity along each LOS. 



C2 Profiles of projected quantities 

From our synthetic maps, we derive radial profiles by aver- 
aging over an azimuthal range of ±45° to both sides of the 
desired direction. The same procedure was used to derive 
the observed profiles (S10). 

Projected profiles are not sensitive to the direction of 
the LOS: e.g. consider maps derived by projection along the 
z-axis and the ?/-axis. For both cases, profiles along the x- 
axis are very similar. 
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C3 Determination of the CF radius 

For a given LOS and direction in a map, we derive the pro- 
jected temperature profile as described in Sect. |C2] We de- 
fine radial ranges with temperature slopes ^ 0.02keV/kpc 
as CFs. This sets an inner and an outer edge for each CF. We 
define the nominal radius of the CF to be centred between 
the front's inner and outer edge. 

The contrast of X-ray brightness or metallicity across 
the CF is the ratio of the this quantity at the inner and 
outer edge of the CF. 



Figure Bl. Time series of temperature slices, taken in the orbital 
plane of the subcluster. White arrows indicate the velocity field 
(w.r.t. main cluster centre), colours code the ICM temperature 
in keV. The series is for the fiducial run M2al00dminl00dmax3 
(high-resolution case). The time is noted in the title of each frame 
along with the scale bar for the velocity vectors. The top left panel 
shows a large-scale view at t = —0.8 Gyr. The remaining panels 
focus on the central 150 kpc. They share the temperature colour 
scale in the bottom left. 
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Figure B2. Profiles of temperature, density, pressure, and metal- 
licity along the x-axis, for the fiducial run (high-resolution case), 
at different timesteps (see legend). For clarity, the pressure pro- 
files of different timesteps are divided by 1.5, 1.5 2 , 1.5 3 , 1.5 4 . For 
the metallicity, we show two sets of profiles (see Sect. |5T) . For 
the upper set, the initial metal density is proportional to the light 
profile of M87. For the lower set, the initial metal density profile 
is flatter. The whole set of flatter profiles is divided by a factor 
of 10 to avoid confusion. On the rhs side (+a>direction), we mark 
the outer CF by a thin vertical doublet line of the colour/style 
matching the legend. 



APPENDIX D: DETAILS OF COMPARISON 

For the sake of brevity, we have omitted some non-essential 
but still useful figures regarding the comparison of simula- 
tions and observations in the main text and collect them 
here. 

Fig. |D1 1 shows the projected profiles of X-ray brightness 
and temperature for the simulation runs that do not match 
the observations. 

Figs. |D2| and |D3| compare brightness residual maps for 
our best-match runs. 
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Figure Dl. Comparison of azimuthally averaged profiles of X- 
ray brightness and projected temperature for runs not matching 
the observation. See Sect. |3. 4] and Fig. [8] for details. 
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Figure D2. Brightness residual maps from the four best- 
matching simulations (first group in Table p] see labels above 
panels). The red circle in the first panel indicates the observed 
field of view. 



This paper has been typeset from a T^X/ WT^X. file prepared 
by the author. 
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Figure D3. Same as Fig. |D2| but for second best simulations 
(second group in Table [51. 



